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I. Introduction 


Entrainment is one of the limits Imposed upon the operation of 
capillary heat pipes. ^ It occurs when the vapor counterflow is fast 
enough to shear off tiny drops of liquid from exposed liquid surface areas 
in the capillary structure. These drops are swept along or entrained with 
the vapor flow. The noisy impingement of these droplets at the condenser 
end of a heat pipe has been cited as the primary evidence of this 
phenomenon. 2 

A modification of the so-called Kelvin-Helmholtz interfacial instabil- 
3 

ity plays a primary role in establishing the conditions for this phenomenon. 
The criterion for avoidance of this instability is often used by designers 
as a conservative estimate of the likelihood of the occurrence of entrain- 
ment in their heat pipes. If 


We 


P V 2 /2 

V V 

y/r o 


< 1 


( 1 ) 


where p v = vapor density, 

= vapor velocity, 
y “ surface tension, 
r Q = capillary pore size, 

entrainment is assumed to be of no importance. One may easily deduce this 

3 

criterion from Lamb's work by making the recognition that the wick structure 
restricts unstable wavelengths to values less than the effective capillary 
pore size. The more common derivation of Eq. (1) is based upon the 
comparison of vapor dynamic pressure and liquid surface tension forces . 
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This familiar background is of interest in connection with the electro- 
hydrodynamic heat pipe because this new device with its extensive free 
surface area is apt to exhibit a modified type of entrainment. The liquid- 
vapor interface is maintained in equilibrium by the electric field 
polarization force; this force thus becomes important in counteracting the 

vapor shear. A simple dimensional analysis might suggest the following as 

4 

a criterion for the avoidance of this instability. 


We 


P V 
v v 


(e.-E )E 
t v o 


< 1 


( 2 ) 


where E q is the tangential electric field intensity measured at the interfaces 
and and are the liquid and vapor dielectric constants, respectively. 

This simple result is , however, not complete, as analysis of a relevant 
electric field coupled interfacial stability problem shows. 
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II. Theory 

A. The Model 

Following methods similar to those of Melcher,^ the general problem 
posed In the Figure 1 is considered. The dielectric liquid (below) 
moves with a uniform velocity profile in the + x direction; the vapor 
(above) moves in the opposite direction, also with a uniform velocity 
profile. This invlscid model, originally attributed to Rayleigh, 
results in a non-allowable flow discontinuity between the liquid and 
vapor flows; large shear forces are present at this Interface which in 
general are not negligible. Based upon this argument significant objection 
may be raised concerning the validity of the model. ^ Still, for the 
purpose of deriving a general design criterion similar to Eq. (1) , the 
Inadequacies of the model are overlooked here. 

The interface is assumed to reside in the region of an electric field, 
mutually perpendicular to the R and y directions. This field is E q at 
the Interface (y => 0) , but exhibits a mild gradient, decreasing in the 
+ ft direction. Mote that purely x-dependent surface perturbations £, 
of the assumed form 

5(x,t) * Re Ue J(wt-kx) ] , |ij « 2 tt/1c (3) 

do not significantly disturb the electric field. The surface is in 
equilibrium at y = 0. 

B. EHD- Coupled Kelvin-He Imho 1 1 z Instability 

The fluid equations for both the liquid and vapor are of similar 


form. 



4 



tangential electric field perpendicular to 
vapor shear flow. 
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- 0 (5) 

The incompressibility assumption for the vapor is valid » due to the low 
vapor velocities considered. Combination of Eq.’s (4) and (5) yields 

V 2 Pl ' = 0 <6) 

for the liquid or vapor perturbation pressure. Consistent with Eq. (3), 
assume 


P^U.y.t) = Re [p 1 (y)e J( “ t_k * > ] (7) 

as typical of the spatial dependence of all perturbation mechanical 
variables. Use of the boundary conditions 


liquid; 

v (x»y * 5_»t) * Re [j(«-kV )£e jC(1)t kx) ] 

y ^ 

(8) 

vapor : 

v y (x ( y = £ + »t) = Re [j (w+kV v Ke^ wt-k *^ ] 

(9) 


yields the result for the perturbation pressure quantities measured at 
the interface. 

P^Cy = S»t) 

P v (y - 5,t) 


k 


5 *“ P A SC 


( 10 ) 


p (oo+kV )‘ 

V V 

k 


C - P v gC 


( 11 ) 


i 
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The surface equation balances the pressure forces against surface 
tension and the electric polarization force. 


P* - S v - k 2 YC -5T e (b 


( 12 ) 


where <5 T (£) is the perturbation electric stress 


S T 


(t r'v )E o 2 l 


<E q > 0) , 


<VS> E o If 


(13) 


(14) 


and the quantity -r=- is always negative for our case. 

% 

The dispersion relation relates the eigenvalues u> and k. It is 
obtained by substituting Eq. (14) into Eq. (12) and canceling out 
Expressed as a polynomial in w, it is 


where 


and 


to + b(k)a) + c(k) = 0 S 


b (k) = - 2k 


" p v V v 


P i + P v 


(15) 


(16) 


c(k) = 


2 2 
k 2 (p t - 


V # + p V ) /p D - P \ 

.A v y . _ f _i ek 

p i + p v W + P v/ 

ic__ k ( V~ .M E |« 

+ p v Vl + p v/ ° 135 


(17) 


Absolute instability occurs if, for some real value or values of k, the 
frequency id exhibits a negative imaginary component. For this case, the 
initially small displacement £(x,t) grows exponentially in time, until 
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some nonlinear mechanism alters or inhibits further growth. Thus, this 
absolute instability criterion does not necessarily predict the onset of 
entrainment, but only the initial amplification of small perturbations from 
which entrainment may result. Such a criterion may be safely used as a 
conservative design factor. 


C. 


The Stability Criterion 

It is helpful to define a new dimensionless quantity. 


We* 





«p r 


P v )g 


<e * " 


e )E 
v o 


JE 

35 


1/2 


(18) 


where V = V + V. . 

v i 

The magnitude of We* is easily shown to determine stability . 


(i) 

(ii) 


He' < 1, stable. 

We* = 1, marginal stability of single wavenumber, 


p * p v V 

k * = 2 y<p; + p v ')" 


(iii) 


He’ > 1, unstable band of wavenumbers. 


(19) 


D. Discussion 

The surface tension continues to play a role in the entrainment limit 
of electrohydrodynamic heat pipes, as evidenced by Eqs . (18) and (19), and the 
stability criterion itself. This conclusion is not surprising, because 
of the nature of the instability itself. As opposed to the nature of a 
capillary wick, the surface of the liquid in an EHD flow structure is 
open, with no restrictions placed upon the wavenumbers of likely 
perturbations. Thus, the wavenumber with the maximum temporal growth rate 
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plays the crucial role in instability. Both surface tension, which 
strongly stabilizes the larger wavenumbers, and the electric field, which 
strongly stabilizes the smaller wavenumbers, must be used to calculate k*. 

If it is found experimentally that entrainment in EHD heat pipes 
presents an unusually severe operating constraint, the means exist to 
alleviate the problem. See Fig. 2. Here the free surfaces along at 
least part of the length of the EHD flow structure of a heat pipe are 
baffled with some insulating wick material (e.g., Refrasil) . 
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vapor 



dielectric interface 
subject to vapor 
shear flow 


vapor 



* ' * * 

dielectric liquid 


“ protective capillary 
baffling 

vapor 


(b) 


Figure 2. Parallel-plate EHD flow structures: 

(a) with dielectric interface subject 
to vapor shear and possible Kelvin- 
Helmholtz instability; (b) with 
capillary baffling. 
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III ♦ Experiment 


A* Experimental Geometry 

The experimental apparatus shown in Figure 3a and b was used to test the 
simple theory described above* A rectangular channel 6" long* 1" wide and 
1/2" deep was equipped with an array of 1/16" diameter brass rods which 
served as electrodes* These rods were flush with the top of the channel and 
were connected with alternate polarity to a high voltage a-c power supply 
(400 Hz). With the liquid surface flush with the tops of the electrode 
rods* an essentially tangential electric field with a negative gradient away 
from the liquid could be imposed. 

The entire channel structure shown in Figure 3 was placed inside a long 
plexiglas enclosure with a variable speed fan at one end. The speed of the 
air flow was controlled by adjusting a rheostat, and it was measured with a 
pitot tube and an inclined manometer. The accuracv of air speed measurements 
was estimated at 10%. 

Care was taken to keep the liquid surface flush with the bottom of the 
wind tunnel. This step minimized the incidence of turbulent eddies and made 
the onset of the Kelvin- Helmholtz instability a more distinct and observable 
transition* See Figure 4. 

B. Procedure 

After filling the channel so that the liquid surface just covered the 
electrode rods (see Fig* 3b), the voltage was turned up to a relatively high 
value. Then the fan was turned on to a predetermined speed. With the voltage 
high enough, the liquid surface remained quiescent, stabilized by the electric 
field gradient force. The voltage was slowly turned down until small amplitude 
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(a) Top view of fluid channel apparatus. 
(Not drawn to scale). 



(b) Cross-sectional view of fluid channel. 

Note that fluid fills reservoir level with 
the top of the channel and with the electrode 
rods . 


Figure 3. Experimental fluid channel used to study 
EHD— coupled Kelvin-Helmholtz instability • 



to manometer 



Figure 4. Wind tunnel used in experiments, 
showing fluid channel placement. 
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wave-like disturbances were seen to develop. This value of voltage, 
measured by an electrostatic voltmeter, was taken as the marginal stability 
point for the given air speed. The voltage was quite distinct and reproducible; 
further decrease of the voltage always produced large nonlinear disturbances 
which washed downstream and spilled out of the channel. Increasing the 
voltage once again always stabilized the surface. 

C. Results 


The data are plotted in Figures 5a and b, which show the relation between 
voltage and air speed at the transition between stable and unstable regimes. 
The theoretical curves shown result from setting We' of Eq. (18) to unity. 
Values above the curve are unstable and values below are stable. Note the 


rather good agreement of the experimental data with the theoretical curve for 
both corn oil and transformer oil.* (The properties of these liquids are given 
in Appendix A.) The small difference between the theory and experiment can be 
accounted for partially by measurement uncertainties, but also by the limita- 


tions of the theoretical model chosen for the Kelvin- Helmholtz instability. 


Further, the value of 


chosen for the theoretical calculation is somewhat 


suspect, due to transverse variation of this quantity at the liquid surface. 


* Due to the tendency of Freon-113 to evaporate rapidly when the air flow 
was introduced, no data for this important dielectric fluid was obtained. 
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Figure 5. (a) Experimental marginal stability data 

with theoretical plot. Fluid is Mazola 
corn oil. 


Air Velocity (m/sec) 
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Figure 5. (b) Experimental marginal stability data 

with theoretical plot. Fluid is G.E. 10-C 
transformer oil. 
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IV. Conclusion 

The success of the theoretical analysis in predicting onset of the Kelvin- 
Helmholtz instability is the best justification of the work reported here. 

Still , the lack of data on freon or other much less viscous dielectric liquids, 
and the relative simplicity of the electrode structure used in these experiments 
Introduce some uncertainty as to the general applicability of the criterion 
developed and tested here to typical EHD heat pipes. The difficulty of 
performing well-controlled laboratory experiments with the flow structures typical 
of an EHD heat pipe preclude a more complete test of the theory at this time. 
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APPEND XX A 

PROPERTIES OF SOME DIELECTRIC FLUIDS 
(0 72° F) 


fluid 

density 

(kg/m 3 ) 

dynamic 

viscosity 

(H - * 2 ) 

\ sec / 

relative 

dielectric 

constant 

surface 

tension 

(N/m) 

Mazola Com 
Oil 

0. 91*10"* 

5.46 *10 -2 

3.1 

0.062 

G.E. 10-C 

transformer 

oil 

0.87*10 3 

1 

~1Q -2 

2.56 

0.04 

Dupont 

Freon-113 

1.57*10 3 

0.68 *10~ 3 

! 

2.41 

0.019 
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APPENDIX B 
Nomenclature 


E 

E 

o 

g 

3 

k 

n 

P 

r 

o 

t 

V 

We 

We E , We' 
y, z 


electric field vector 

tangential electric field at fluid surface 

2 

gravitational acceleration (9.81 m/sec ) 
wavenumber 

unit vector normal to fluid surface 
fluid pressure 

effective capillary pore size 
time 

fluid velocity 
capillary Weber number 
electric Weber numbers 
coordinate axes 


Y 

£ 

0 

U) 

5 


surface tension 
dielectric constant 
fluid density 
radian frequency 
fluid surface displacement 


% — liquid 


Subscripts : 


v 


vapor 
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